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ABSTRACT 

We study the suprathermal electron acceleration mechanism in a perpendicu- 
lar magnetosonic shock wave in a high Mach number regime by using a particle- 
in-cell simulation. We find that shock surfing/surf atron acceleration producing 
the suprathermal electrons occurs in the shock transition region where a series 
of large amplitude electrostatic solitary waves (ESWs) are excited by Buneman 
instability under the interaction between the reflected ions and the incoming 
electrons. It is shown that the electrons are likely to be trapped by ESWs, 
and during the trapping phase they can be effectively accelerated by the shock 
motional/convection electric field. We discuss that suprathermal electrons can 
be accelerated up to mjC 2 (t>o/c), where m^c 2 is the ion rest mass energy and vq 
is the shock upstream flow velocity. Furthermore, some of these suprathermal 
electrons may be effectively trapped for infinitely long time when Alfven Mach 
number Ma exceeds several 10, and they are accelerated up to the shock potential 
energy determined by the global shock size. 

Subject headings: acceleration of particles - cosmic rays - plasmas - shock waves 



1. Introduction 

The origin of high energy particles seen radiating in astrophysical synchrotron sources is 
still a major unresolved problem of high-energy and plasma astrophysics. These sources in- 
clude supernova remnants (SNRs) and extragalactic radio sources by jets etc. (e.g., Koyama 
et al. 1995; Meisenheimer et al. 1989; Carilli and Barthel 1996; Junor et al. 1999). The 
shock acceleration has been discussed as one of important processes producing the high- 
energy particles, and there are many theoretical and observational efforts on the high-energy 
particle acceleration/heating so far. The diffusive /Fermi acceleration process particularly 
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met with great success (Axford et al. 1977; Bell 1978; Blandford and Ostriker 1978), because 
it predicts a power-law energy spectrum with the power-law index of 2 that depends weakly 
on the shock compression ratio, and because the index suggests a very close to that believed 
to exist in many astrophysical sources. 

From in situ measurements at the planetary shocks and at the interplanetary shocks, 
the energetic ions with a power-law spectrum has been observed (e.g., Scholer et al. 1980; 
Gosling et al. 1981), and they are believed to be produced by the diffusive shock acceleration. 
For the energetic electrons in the astrophysics context, the radio spectra of SNRs and the 
emission from active galactic outflow, for example, are believed to be produced by the shock 
acceleration, and much work has been concerned with the diffusive shock acceleration of 
electrons. However, the energetic electrons with a power-law spectrum at the interplanetary 
shocks and at the planetary bow shocks are rarely ever observed (e.g., Shimada et al. 1999), 
and the diffusive-type acceleration seems to be malfunctioning in the Heliosphere. The 
electron nonthermal acceleration still remains an unresolved issue of considerable interest. 

In order to the shock diffusive acceleration to start, a certain pre-acceleration mechanism 
should provide a seed population having from thermal to mildly high energy, and then 
the energized particles can be subsequently accelerated to highly relativistic energies. Our 
understanding of the electron pre-acceleration mechanism, which is the so-called injection 
process, is very limited (e.g., Levinson 1996). Furthermore, MHD waves which provide the 
resonant scattering of electrons in the shock upstream are not well studied so far. 

By taking into account plasma instabilities in detail, Papadopoulos (1988); Cargill 
and Papadopoulos (1988) proposed the electron energization process at high Mach num- 
ber shocks, in which electron heating is produced through two-step instabilities in the shock 
transition resion where the reflected ions coexist with the incident ions and electrons (Leroy 
et al. 1982; Wu et al. 1984). Buneman instability (BI) (Buneman 1958) is first excited by the 
velocity difference between the reflected ions and the incident electrons, and the electrons 
are heated up by the instability. As the next step, the ion acoustic (IA) instability is trig- 
gered under the pre-heated electron plasma by BI, and the electron temperature increases 
by a factor approximately equal to M\j ' (3 e , where M A is Alfven Mach number and f3 e is 
the upstream electron plasma beta. This scenario of the electron heating mechanism has 
been demonstrated by using a hybrid simulation code where ions are treated as particle but 
electrons are assumed to be massless fluid (Cargill and Papadopoulos 1988). Recently Shi- 
mada and Hoshino (2000) extended their study to the shock structure including the electron 
dynamics by using a particle-in-cell simulation where both ions and electron are treated as 
particle. They found the formation of the phase space "electron hole" in association with 
a localized, large-amplitude, electrostatic solitary wave in the shock front region. It is dis- 
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cussed that BI in the nonlinear stage produces a large-amplitude, localized electrostatic wave 
which is called as the electrostatic solitary wave (ESW) (Davidson et al. 1970; Omura et al. 
1994). They also found that the ESW itself is almost stable, but the interaction with other 
waves in the inhomogeneous shock transition region leads to the dynamical evolution of the 
shock energy dissipation. They discussed that the electrostatic solitary wave (ESW) plays 
an important role not only on the rapid electron thermalization but also on the nonthermal 
electron acceleration, and that the ESWs obtained by the simulation can be compared with 
in situ observations of the electric field waveform at the Earth's bow shock (Matsumoto et al. 
1997; Bale et al. 1998). However, the actual physics of the nonthermal electron acceleration 
processes has remained elusive. 

In this paper, we study in details the suprathermal electron acceleration through the 
interaction of electrons with ESWs in a transverse, magnetosonic, high Mach number shock. 
We discuss that the electrons trapped by ESW can resonate with the shock motional electric 
field, and the so-called shock surfing mechanism is effective for producing the non-thermal, 
high-energy electrons (Hoshino 2001; McClements et al. 2001). Shock surfing acceleration at 
quasi-perpendicular shocks is usually considered for ions to be a pre-acceleration mechanism 
to initiate diffusive shock acceleration (Sagdeev 1966; Sagdeev and Shapiro 1973; Sugihara 
and Mizuno 1979; Katsouleas and Dawson 1983; Lembege and Dawson 1984; Ohsawa 1985; 
Zank et al. 1996; Lee et al. 1996; Ucer and Shapiro 2001). The electrostatic field at the 
shock front directs upstream, and the ions which energy cannot overcome the electrostatic 
potential well are reflected upstream in a super-critical shock of M A > 2.7. In the standard 
shock surfing, the ions are trapped between the shock front and the upstream by the Lorentz 
force. During the reflection process, ions travel along the shock front and can be accelerated 
by the motional/convection electric field. In this mechanism, however, electrons can be 
neither reflected nor accelerated. We discuss that a series of large amplitude electrostatic 
waves excited by the nonlinear Buneman instability in the shock transition region can play 
a role of counterpart of the shock front potential for the case of the ion shock surfing, and 
those ESWs can effectively trap electrons. Then the electron shock surfing mechanism can 
be switched on. 

We also propose that for a large Mach number shock with Ma > several 10, electrons are 
likely to be trapped for infinitely long times, because the electrostatic force becomes always 
larger than the Lorentz force. In this situation, the electrons can be perfectly trapped by 
ESWs, and can be quickly gain their energies until they obtain the global transverse shock 
potential energy. We discuss that the shock surfing/ surfatron near the shock front region is 
the efficient accelerator for the relativistic electrons. 
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2. Nonthermal Electrons in High Mach Number Shocks 

We study the ion-electron dynamics organizing the electron energization and wave activ- 
ities in the shock transition region by using the one- dimensional, particle-in-cell simulation 
code where both ions and electrons are treated as particle (Hoshino et al. 1992). In our sim- 
ulation study, a low-entropy, high-speed plasma consisting of electrons and ions is injected 
from the left boundary region which travels towards positive x. At the injection boundary 
at x = 0, the plasma carries a uniform magnetic field B z , polarized transverse to the flow. 
The downstream right boundary condition is a wall where particles and waves are reflected. 
Under the interaction between the plasma traveling towards positive x and the reflected 
particles and waves from the right-hand boundary, the shock wave is produced, and it prop- 
agates backward in the —x direction. Then the downstream bulk plasma speed becomes 
zero in the simulation frame. The grid size is comparable to the electron Debye length in 
downstream, includes about 1800 particles for each species in downstream. The plasma 
parameters are as follows; upstream plasma f3 e — fa — 0.01 (f3j = 8imTj/B 2 , where n, T, 
and B are, respectively, the density, temperature, and magnetic field strength), the ratio of 
the plasma frequency to the electron cyclotron frequency u pe /uj ce = 19 (uj pe = a/ 'Arme 2 / 'm e , 
u ce = eB/(m e c)), and the ratio of ion to electron mass m,i/m e = 20. The whole system size 
L/(c/uj pe ) = 300, the upstream speed vq in the simulation frame is 0.25c, and the ratio of 
Alfven speed to the speed of light is Va/c = 1.2 x 10~ 2 . The shock Alfven Mach number 
becomes Ma ~ 32 in the shock frame. (Note that the upstream speed in the shock frame is 
about (3/8)c.) Our simulation is performed until tui pe = 1815, when the shock front reached 
at X/(c/uo pe ) ~ 100. 

Figure 1 shows a snapshot of the nonlinear stage of a collisionless shock front region at 
tujp e = 735. The left-hand panel shows the shock front region from X/(c/u pe ) = 150 to 300, 
while the right-hand panel is its enlarged picture. The leftward and the rightward regions 
are respectively the shock upstream and downstream. From the top, the ion phase space 
diagram in (X, U ix ), the electron phase space diagrams in (X, U ex ), and in (X, \U e \) where 
\U e \ = (U 2 X + U 2 y + U 2 Z ) X I 2 . The bottom three panels are the transverse electric field E y , 
the magnetic field B z , and the electric field E x . The plasma four velocity U is normalized 
by the upstream flow velocity U$ = Vq/ 'y/l — (v /c) 2 . The magnetic field, and the electric 
field are normalized by the upstream magnetic field Bq and the upstream motional electric 
field Eq = vqBq/c, respectively. The spatial scale is normalized by the electron inertia length 
c/u pe in upstream. 

In the shock front region in Figure 1, we can find two ion components in the ion phase 
space diagram; one is the cold ion flowing into the shock downstream, the other is the 
reflected ion going away from the shock front. This ion dynamics has been well studied by 
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both the satellite observations (Paschmann et al. 1981; Gosling et al. 1982; Sckopke et al. 
1983) and the theory/simulation studies for super-critical Mach number shocks (Leroy et al. 
1982). A part of the incoming ions is reflected due to both the polarization electrostatic field 
and the compressed magnetic field in shock downstream. The magnetic overshoot structure 
can be found around X/ (c/u pe ) ~ 248. 

In addition to such ion dynamics around this shock transition region, we find the elec- 
tron hole in the electron phase space diagram in (X, U ex ) around X/(c/uj pe ) ~ 223, and its 
corresponding electrostatic solitary wave (ESW) in the electric field E x . Other ESWs in the 
growing phase can be also found in the foreside of the largest ESW structure. Shimada and 
Hoshino (2000) discussed that those localized structures are excited by Buneman instability 
between the reflected ions and the incoming electrons (Papadopoulos 1988; Cargill and Pa- 
padopoulos 1988; Dieckmann et al. 2000a,b). The simulation result in Figure 1 is basically 
same as that obtained by Shimada and Hoshino (2000) and Schmitz et al. (2001). 

Looking at the electron hole structure in the shock transition region in details, we 
find that some electrons are accelerated and gain a large amount of energy. After the 
strong acceleration in the shock transition region, the electrons are transported downstream. 
Figure 2 shows the downstream electron energy spectrum. The energy 7 = a/1 + (U e /c) 2 is 
normalized by the incident electron bulk energy 70 = 1/ a/1 — (v /c) 2 ~ 1.03. The spectrum 
is superposed over 538 snapshots for the time interval from u pe t = 535 to 1815. The dot- 
dashed line is the Maxwellian fit for the spectrum, and the bottom dotted line is the so- 
called one-count level, which is the reciprocal of the number of snapshots of 1/538 in this 
case. Although the statistics of the nonthermal high energy tail is not good, but we can 
clearly find the enhancement of the nonthermal population above the Maxwellian level over 
7/70 > 1.8. 

Shimada and Hoshino (2000) also obtained that the similar suprathermal energy spec- 
trum for the high Mach number shock of Ma ~ 10. In addition, they found that the 
suprathermal electrons are not generated for a low Mach number case of Ma ~ 3, and the 
spectrum is well fitted by a thermal Maxwellian. With regard to the electron energy spec- 
trum obtained by these simulations, a high Mach number shock seems to be a candidate for 
producing high energy electrons. 

Before discussing the mechanism producing the suprathermal electrons, we would like 
to give a comment that the suprathermal electron energy is much larger than the potential en- 
ergy of ESW. The potential energy ec/> esw = eE esw A esw can be estimated as 20(u ce /u pe ) (v /c)m e c 2 ~ 

0. 26m e c 2 , where we have assumed E esw ~ 20E , the width of ESW, A 

csw ~ c/ ui pe from Figure 

1. Other simulation parameters are u ce /uj pe = 1/19 and vq/c = 1/4. Then we get the incre- 
ment of the Lorentz factor normalized by the initial Lorentz factor as Aj/^q ~ (7 — 7o)/7o ~ 
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0.26. Therefore the electrons cannot be accelerated over 7/70 ~ 1.26. Moreover, due to a 
bipolar electric field structure of ESW, no energy gain is expected after crossing the whole 
ESW in the x-direction. 



3. Electron "Shock Surfing Acceleration" 

In order to know how and where the high energy electron is accelerated in a high Mach 
number shock, we study the electron trajectory near the shock front region. Figure 3 shows 
the stack plot of the E x field and two typical electron orbits (left-hand panel) and the time 
history of the total momentum \P e \ — m e \U e \ (right-hand panel). The total momentum \P e \ 
is normalized by the upstream momentum P = mJJ^. The vertical axis t is normalized by 
the electron plasma frequency u pe , and the time history is plotted after toUJ pe = 535. The 
most of waveforms in the shock transition region are basically propagating the negative x 
direction, which speeds are slightly faster then the shock front propagation speed in the wave 
growing phase, because ESWs are generated by the velocity difference between the reflected 
ions and the incoming electrons, and because the reflected ions carry larger momentum than 
the incoming electrons. As the wave amplitude increases in the nonlinear phase, the speed 
of ESW becomes slow down, and the waveform is convecting downstream. We find that 
some waves are sporadically transported downstream during (t — t )uj pe = 200 ~ 420. Note 
that the downstream plasma is at rest in our simulation, and the shock front is propagating 
towards the negative x direction. 

Let us first look at the trajectory denoted by the solid line (left-hand panel). The 
particle is situated in the upstream region before (t — to)ou pe ~ 100, and it is convected 
towards the shock front with the E x B convection speed of Vq = cE /B . When the particle 
arrived at the shock front at (t — t )u pe ~ 120, it is trapped by a large amplitude ESW 
until (t — t )uj pe ~ 220, and after that it escapes out of the ESW structure and is convecting 
downstream. From the time history diagram (right-hand panel), we find that the efficient 
energization occurs while the electron is trapped for the time interval from (t — t )u pe = 120 
to 240. The gradual energy increase can be also found after the detrapped phase from 
(t — to)u> pe = 240 to 500. The second phase of acceleration corresponds to the region having 
the magnetic field gradient in association with the magnetic overshoot. 

The dashed line shows another typical electron orbit for which the interaction of electron 
with ESW was very weak. When crossing the shock front region, the electron was not 
trapped by ESW, and it passed through ESW without any significant energy gain. The 
sinusoidal oscillation means the gyro-motion of particles. Note that the solid curve has a 
larger sinusoidal oscillation than the dashed curve, because the electron denoted by the solid 
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curve gains larger energy during the shock front crossing. 

Figure 4 shows the time history of the first adiabatic invariant of P ±e / B 2 for the same 
particles analyzed in Figure 3. (P± e /B 2 ) is averaged over the upstream electron gyro-period 
of cj" 1 , and is also normalized by the upstream value of P± /Bq. The dashed curve is 
almost constant, and we know that the electron motion is almost adiabatic during the shock 
crossing. For the solid curve case that had the strong interaction of electron with ESW, 
we find a sharp energy increase for the time interval from (t — to Vpe — 120 to 240, which 
is suggestive a non-adiabatic acceleration. After (t — to)u pe ~ 240 we found the gradual 
energy increase in Figure 3b, but the baseline of the normalized invariant seems to be almost 
constant except for a sinusoidal oscillation. This suggests that the electron heating is almost 
adiabatic after the electrons are convected downstream, and the non-adiabatic acceleration 
occurs only when the electron interacts with ESW in the shock transition region. 

Figure 5 shows the above two particle trajectories in the x — y plane. Both two trajec- 
tories are drawn from (X, Y) ~ (200 c/u pe , 0) at u pe t = 535 in the upstream region, and the 
electrons are convected towards positive x along y — 0. The Larmor radii are very small in 
the upstream region because their thermal velocities are very small. Note that this result 
was obtained by the one-dimensional, particle-in-cell simulation, but both positions of x and 
y are calculated. Also note that the x and y distances are plotted on the different scale in 
Figure 5. They encounter the shock front region around X/(c/uj pe ) = 225. The notation of 
the solid and dashed curve is same as that in Figures 3 and 4. We find that the electron 
denoted by the solid curve are traveling along the y axis in the shock front region, it can be 
understood that the main energy gain of the electrons comes from the motional electric field. 
In fact, the energy gain of electron from the motional electric field Ae m can be estimated 
as eE m AY, where AY ~ Q0(c/uj pe ) and E m = (v — v esw )B /c is the motional electric field 
in the ESW frame. Since the propagation speed of ESW is almost same as the speed of 
the reflected ions (Shimada and Hoshino 2000), we may assume f esw ~ — ^o- Then we get 
Ae rn ~ \1Q{ui ce / ui pe ) {vq/ c)Tn e c 2 ~ 1.6m e c 2 , and the Lorentz factor of the accelerated electron 
becomes 2.6. Namely, the normalized momentum P e /Po ~ 9, which is almost consistent with 
the momentum gain in Figure 3. 

From the above analysis of the simulation data, we think the so called "shock surfing" 
mechanism plays an important role on electron acceleration. Let us quickly review the idea of 
the shock surfing (Sagdeev 1966; Sagdeev and Shapiro 1973; Sugihara and Mizuno 1979). The 
shock surfing mechanism has been extensively studied for the ion acceleration (Katsouleas 
and Dawson 1983; Lembege and Dawson 1984; Ohsawa 1985; Zanket al. 1996; Lee et al. 1996; 
Ucer and Shapiro 2001). Due to the inertia difference between ions and electrons flowing into 
the shock, the polarization electric field normal to the shock front is formed. An ion having 
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a small velocity normal to the shock front, i.e., l/2Mvl < e0 s , can be reflected from the 
shock front to the shock upstream, where <p s is the electrostatic shock front potential induced 
by the inertia difference between ions and electrons. During the gyro-motion in upstream, 
ions gain their energy from the motional electric field parallel to the shock front. If the 
width of the shock front potential S is thin, the electric force can overcome the Lorentz 
force, and then the multiple reflection can occurs. In a quasi-perpendicular shock, Ohsawa 
and Sakai (1987) studied the nonlinear steepening of magnetosonic waves, and showed that 
if u pe /u ce < 1, the thickness of the shock front potential (f> s becomes of the order of the 
electron inertia scale c/u pe . They then demonstrated by using a particle simulation the 
strong resonant interaction of the reflected ions with the motional electric fields. 

What we discuss here is the "electron" shock surfing acceleration (Hoshino 2001; Mc- 
Clements et al. 2001). The above "standard" shock surfing acceleration cannot apply for 
the electron acceleration, because the electron cannot be reflected from the shock front by 
the shock front potential <p s . We propose a new scheme of "electron" shock surfing accel- 
eration under the action of ESW. Since the electron hole is a positively charged structure, 
and an electron can be trapped if m e v x /2 < e<p esw , where e sw is the scalar potential for the 
electrostatic solitary wave (ESW). (Note that the Lorentz force is also an important agent 
for the particle trapping, but we will discuss this point later.) Furthermore, the propagation 
velocity of ESW differs from the plasma bulk velocity, and ESW together with the trapped 
electrons can stay longer time in the shock transition region. Namely, in the frame moving 
with ESW, the convection electric field is not zero. Therefore, we think that the so called 
"shock surfing" mechanism is occurring for electrons. 

Figure 6 summarizes our idea of the electron shock surfing mechanism. Top panel shows 
an electron's trajectory in the x — y plane. The magnetic field is polarized perpendicular to 
the x — y plane, and the plasmas are convecting towards positive x. The shock upstream 
is the left-hand side, while the downstream is the right-hand side. Bottom panel shows the 
electric field E x along the x axis, and ESW in association with its electron hole in phase space 
is depicted in the center. Due to the nature of the electron hole, the electron charge density 
is slightly lower than the ion one, and ESW has a bipolar signature with diverging electric 
field. If an electron convecting towards the ESW structure is reflected by both the Lorentz 
force and the electric field E x and is trapped inside the ESW structure, it is successively 
accelerated towards the negative E y direction. As increasing the electron's velocity v y by 
the shock surfing/surfatron acceleration, it can be de-trapped from ESW when the Lorentz 
force ev y B z /c becomes larger than the electric force eE x , and then it is convecting towards 
downstream and becomes quickly an isotropic, gyrotropic distribution. The time variable 
nature of the ESW evolution may be the other important factor to control the de-trapping 
process. 
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4. Unlimited Electron Acceleration 

To illustrate the efficiency of the above electron shock surfing mechanism, let us discuss 
the x component of the equation of motion, 

d e 

—P x = eE x + -v y B z , (1) 

where P x = m e U ex is the electron momentum. In our interest, E x on the right-hand side can 
be replaced by the electrostatic solitary wave (ESW) produced in the shock transition region. 
If the electric force of eE esw is larger than the Lorentz force of (e/c)v y Bo, the electrons are 
trapped and gain their energy. During the non-adiabatic acceleration phase, the velocity v y 
increases. If the electron satisfies the condition of eE csw < (e/c)v y B , it can escape from 
ESW and will be convected downstream. 

The amplitude of ESW may be estimated by equating the wave energy density to the 
drift energy density of the incoming electron (Ishihara et al. 1981), and we obtain, 

^ ~ \m e nV^ (2) 

where Vd ~ 2vq is the relative velocity between the reflected ions and the incoming electrons. 
The energy conversion factor a is of order of 0(1). Although the nonlinear saturation process 
of BI still remains controversial, we adapt here the conversion factor a discussed by Ishihara 
et al. (1981); Dieckmann et al. (2000b), and we use, 

« = 7 ~ (— ) 1/3 - (3) 
4 mi 

Before discussing the efficiency of electron acceleration, it is better to check whether 
or not the above estimation is reasonable. We first analyze the magnetospheric observation 
data for a high Mach number shock. From the satellite observations in the Earth's bow 
shock, we know that E csw = 50 ~ 200 mV/m (Matsumoto et al. 1997; Bale et al. 1998), 
while the solar wind, motional electric field E is a few mV/m. Then, the ratio of E esw /E 
is of order of 10 1 ~ 10 2 . On the other hand, by virtue of Eq.(2), the ratio of E esw to the 
motional electric field E = v B /c in upstream can be obtained, 



e v A y mi 1 ; 

Since the Alfven velocity Va is about 300 km/s in the solar wind, the theoretical estimation 
of Ecsw/Eq becomes 13 ~ 23. This value is remarkably close to the observation. In the same 
way as this, we can check this for our simulation result. We observed £ , esw /£ , = 10 ~ 30 in 
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Figure 1, while Eq.(4) suggests E esw /E = 19 ~ 23 for m;/m e = 20 and V A /c = 0.012. The 
comparison between the theoretical estimation and the simulation gives a good agreement 
again. Therefore, we think the estimation of the amplitude of ESW in Eq.(2) is reasonably 
well. 

Now we can estimate the maximum electron energy. The electric force Fe = eE x can 
be given by, 

F E ~ eE a „ ~ 2 eB „^£, (5) 

and the Lorentz force Fb in the shock transition region is, 

F B ~ % y B . (6) 
By equating Fe to Fb, we obtain the maximum velocity for v Ytmax as, 

j 2cM A for 2M Ay /a(m e /mi) < 1 

z for 2M / iy / a(m e /m i ) > 1 



If M A > (l/2) v /m 1 /(m e a) ) then i*£ > .Fg is always satisfied, and the maximum speed 
becomes the speed of light. In this case, the electrons cannot escape from the ESW region, 
and the electrons continue to accelerate until they reach the edge of a global shock structure. 
The condition of the "unlimited electron acceleration" for the real mass ratio of mi/m e = 
1830 is given by, 

M A > 43 ~ 75. (8) 

We think that the electron surfing with ESW in the shock transition region is a strong 
accelerator. 



5. Conclusions and Discussion 

We discussed the electron shock surfing through a large amplitude, small-scale electro- 
static field in the shock transition region. The electrons may be unlimitedly accelerated when 
Alfven Mach number M A > 43 ~ 75. It is well known that the dynamical shock emerges 
for a super-critical shock with Alfven Mach number M A > 2.7, and a fraction of incoming 
ions are reflected from the shock front. Due to the velocity difference between the incoming 
electrons and the reflected ions, a series of large amplitude electrostatic waves are excited by 
the nonlinear Buneman instability. We found that the electrons can be trapped by the large 
amplitude wave and the non-adiabatic acceleration of electron occurs under the interaction 
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with the shock motional electric field parallel to the shock front. We also found that, as 
further increasing the Mach number and if M4 exceeds 43 ~ 75, the electrons can be 
unlimitedly accelerated and gain the relativistic energy. Such high Mach number shocks can 
be found in SNRs, which are known as the relativistic electron accelerator with a strong 
X-ray emission. 

It is an open question whether or not the shock surfing acceleration discussed in this 
paper can produce the power-law energy spectrum. For the case of the ion shock surfing, 
Zank et al. (1996) pointed out that shock surfing can produce the non-thermal power-law 
energy spectrum in the shock front region. It may be interesting to study the similar scheme 
for the electron shock surfing. Our understanding on the electron shock surfing is still limited, 
but we would like to emphasize that the acceleration rate of the shock surfing is much faster 
than that of the shock diffusive /Fermi acceleration, and it may be important not only for 
the injection process of the shock diffusive acceleration but also for high energy particle 
acceleration. 

The unlimited electron shock surfing acceleration is based on the assumption that ESW 
stably exists during the interaction. It is known, however, that the shock front is not stable 
and periodically breaks. Namely, the shock reforms with a cycle time of ion gyro-poeriod, 
u C i, which is called as the shock reformation process (Quest 1986). The reflected ion, which 
can excite ESWs through Buneman instability, has also the same dynamical behavior with 
the shock reformation, and BI may be quenched every shock reformation cycle. Therefore, 
most of electrons may be de-trapped from ESW in association with the shock reformation 
process. Assuming that the ion reformation time scale is about ion gyro-period, u;" 1 , the 
maximum energy of the electron shock surfing might be limited by, 

eE c 2 v Q 
£max ~ ~ mc {—). (9) 

0J c i c 

In the case of the simulation parameter used in Figure 2, the above maximum energy nor- 
malized by the upstream flow energy can be estimated as 7max/7i ~ (n r ii/n r ie)(vo/c) ~ 7.5. 
However, the electron energy density became below the one-count level around 7/70 ~ 4.5 
due to the limited number of particles in the simulation box. It seems that we have to 
increase the number of simulation particle at least 100 times more to examine the cut-off 
signature, and we need the larger computer resources to understand the behavior of high 
energy spectrum. 

Our simulation was performed in one-dimensional system, but the multi-dimensional 
effect may be important for the non-thermal particle acceleration. One may think that ESW 
has a finite transverse scale in two- or three-dimensional system, and that the finite ESW 
structure may limit the electron shock surfing acceleration. On the other hand, it may be 
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possible that some electrons move from an ESW to another ESW in a stochastic manner and 
resonate successively with the motional electric field, because the shock surface is expected 
to have a wavy structure. In this case, the maximum energy of the shock surfing is not 
limited by Eq.(9), and we think that the electrons can be unlimitedly accelerated until they 
reach the edge of a global shock structure. We are planning to study two-dimensional shock 
dynamics as the future problem. 

The propagation speed of the electrostatic solitary wave against the shock front is also 
an important agent to control the electron surfing acceleration. If ESW is moving with the 
same speed of the plasma medium, the shock motional electric field E m in the ESW frame 
disappears, and no surfing acceleration occurs. In fact, we find that the efficient acceleration 
occurs in the upstream side of the shock transition region where the amplitude of ESW 
is not necessarily large. ESWs are growing their amplitudes towards downstream in the 
transition region, but the ESW structure starts to move with the almost same velocity with 
the plasma medium, and the motional electric field in the ESW frame becomes small. Then, 
the electrons cease the surfing acceleration. However, ESWs start to collapse by colliding 
with other waves in inhomogeneous shock transition region (Shimada and Hoshino 2000), and 
those collapsed waves may provide further stochastic and surfatron acceleration, as discussed 
recently by McClements et al. (2001). 

In this report, we restrict our attention to the weakly magnetized shock where the 
electron cyclotron frequency u ce is smaller than the electron plasma frequency u pe , because 
most astrophysical plasmas have a weakly magnetized plasma with uj ce < u, pe . In the case 
of u ce ^> uj pe , Buneman instability and its subsequent nonlinear evolution discussed in this 
paper may be suppressed (Dieckmann et al. 2000a). Bessho and Ohsawa (1999) studied 
a quasi-perpendicular magnetosonic shock with u ce ^> u pe using a particle-in-cell simula- 
tion code, and investigated acceleration of electrons during the nonlinear decay process of 
the magnetosonic wave. It is interesting to know the efficiency of suprathermal electron 
acceleration as the function of uj ce /uj pe . 

Finally, we should like to make the point that it may be possible to explore many aspects 
of the electron shock surfing acceleration from in situ observations at the interplanetary 
shocks and at the planetary bow shocks. It is known that the interplanetary shocks initiated 
by the solar flare may develop into a high Mach number shock with Ma ~ 40 around 
the Mercury orbit of ~ 0.4 AU (e.g., Smart and Shea 1985). Such observations could 
bridge a gap of our understanding between the electron injection and acceleration at the 
interplanetary/planetary bow shocks with relatively low Mach numbers and the astrophysical 
shocks with very high Mach numbers. 
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Fig. 1. — High Mach number, magnetosonic shock structure for Ma = 32: an overall shock 
front region (left), and an enlarged picture in the shock transition region (right). 
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Fig. 2. — Downstream suprathermal electron energy spectrum. Dotted lines indicate 
Maxwellian distribution. 
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Fig. 3. — Time evolution of wave form for the longitudinal electric field E x and two typical 
particle trajectories (left), the time history of the total momentum (right). 




Fig. 4. — The time history of the 1st adiabatic invariant P e j 1 B 2 for the two typical particles 
shown in Figure 3. 
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Fig. 5. — The particles trajectories in the plane perpendicular to the magnetic field B z for 
the two typical particles shown in Figure 3. 
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Fig. 6. — An illustration of the electron shock surfing mechanism under an action of ESW 
(electrostatic solitary wave). 



